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The status of the Monte Carlo programs for the simulation of the r-lepton production and decay in high energy 
accelerator experiments is reviewed. In particular, the status of the following packages is discussed: (i) TAUOLA for 
r-lepton decay and PHOTOS for radiative corrections in decays, (ii) MC-TESTER for universal tests of the Monte 
Carlo programs describing particle decays, (iii) KORALB, KORALZ, KKMC packages for r-pair production in e + e~ 
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1. INTRODUCTION 



The package TAUOLA |jg,|3j,y for the simulation 
of rdepton decays and PHOTOS J|,|| for the sim- 
ulation of radiative corrections in decays, have 
a rather long history. Written and maintained 
by well defined authors, they nonetheless mi- 
grated into a wide range of applications, where 
they became ingredients of complicated simula- 
tion chains. As a consequence, a large number of 
different versions are currently in use. From the 
algorithmic point of view, they often differ only 
in a few small details, but incorporate substantial 
amounts of specific results from distinct r-lepton 
measurements. Such versions were mainly main- 
tained by the experiments taking precision data 
on r-leptons. On the other hand, many new ap- 
plications were developed recently, often requir- 
ing program different versions because of inter- 
faces to other packages. 

In the following, I will concentrate on those top- 
ics where changes with respect to the status pre- 
sented at the Victoria r conference two years ago 
(7) were introduced. Since that time, there were 
no changes introduced into the PHOTOS Monte 



"This work is partly supported by the Polish State Com- 
mittee for Scientific Research (KBN) grants 2 P03B 001 
22, and also by the European Community's Human Po- 
tential Programme under contract HPRN-CT-2000-00149 
Physics at Colli ders. 



tHome page at http://wasm.home.cern.ch/wasm/ 



Carlo functionality, and also the TAUOLA inter- 
faces to KORALB Hfl, KORALZ @, and KKMC (ll) 
remain unchanged. On the contrary, the univer- 
sal interface of TAUOLA evolved, and new applica- 
tions, in particular observables for the measure- 
ment of the Higgs-boson parity, became possi- 
ble. The new program MC-TESTER |l| instru- 
mental in the development of future versions of 
TAUOLA, was developed. A new choice of hadronic 
currents for r — > Airv decay modes became avail- 
able. It is based on Novosibirsk data. 

Let me concentrate, in the following three sec- 
tions, on these topics, and close my contribution 
with a summary. 

2. THE r-LEPTON FOR HIGGS BOSON 
PARITY 

In many applications, the precision of the con- 
trol of the r spin effects is not crucial. This is 
the case in searches of new particles yet to be dis- 
covered, or in applications where r-lepton decays 
contribute as a final state of some rare decays of 
known particles. Nonetheless, in such cases, spin 
effects can be of some use as well. As these, there 
is no motivation to develop sophisticated spin al- 
gorithms for every individual case: less precise, 
but universal solutions are welcomed. 

The universal solution of p3[ , based on the 
HEPEVT common block of F0RTRAN77 is now dis- 
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tributed with TAUOLA. The basic idea was to 
calculate the spin (helicity) state of the decay- 
ing r from the kinematical information avail- 
able in HEPEVT and some very simple assump- 
tions on the production mechanisms. The pro- 
gram was checking if the production was through 
// — > Z/7 — > r + r~, W or Higgs-boson inter- 
mediate states/processes. The main properties of 
the algorithm were presented in Victoria already 
2 years ago; modifications necessary for the pro- 
gram to work in the case of some r production 
processes involving new, to-be-discovered parti- 
cles, have gradually been introduced, but one can 
sec [[l4| that not in all cases is the effort com- 
pleted. 

In the meantime algorithm functionality was 
also extended: full spin effects were introduced in 
the case of Higgs-boson decay fL5[ . The case of 
the Higgs boson is exceptionally easy, since the 
full density matrix of the r-lepton pair produced 
from a Higgs-boson is fully defined by boson par- 
ity and r leptons four-momenta. 

This algorithm helped us to design an observ- 
able |l6| for the possible measurement of Higgs 
boson parity in future accelerators such as Lin- 
ear Colliders or the LHC. Let us recall the main 
principle of this technique, which from the point 
of view of spin analysis is quite involved, as it 
requires a study of full spin correlation spanning 
over three levels of decay cascade: h/A — > r + r~, 
t ± — > p ± i> and / o ± — > ^it . The presence of 
two non-observable neutrinos as well as our in- 
ability to reconstruct the Higgs-boson rest frame 
sufficiently well (with a precision comparable to 
the r-lepton mass) complicates the picture even 
further. 

We started from the observation that the dis- 
tribution of the acoplanarity angle of two planes 
spanned on decay products of p + — > 7r + 7r° and 
p~ — ► 7r~7T° (defined in the rest frame of a p + p~ 
pair) is quite sensitive, see fig. jjj to the parity 
of Higgs boson. One has to select sub-samples of 
particles, with the same sign (opposite sign) of 
the energy difference between ir + and ir° in the 
r + rest frame and that between ir~ and ir° in the 
t~ rest frame. This condition is unfortunately 
impossible to realize in practice, as r-lepton mo- 
menta cannot be reconstructed. 
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Figure 1. Acoplanarity distribution of the p + p~ 
decay products ' in the rest frame of the p + p~ pair. 
A cut on the 7r ± to tt° energy difference defined, 
in the generator level t^ 1 rest frames, to be of the 
same sign is used. No smearings included. The 
thick line denotes the case of the scalar Higgs- 
boson and the thin line the pseudoscalar one. 

In the next step of our study we have weakened 
the requirement, and we have used a replacement 
for the r-lepton momentum; the p + and p~ di- 
rections (in the rest frame of a p + p~ pair) re- 
place the unobservable directions of r flights and 
the Higgs mass constraint provides the r-lepton 
energy (also in the same frame). Indeed, with 
such a technique, and after including some as- 
sumption on detector smearings, the sensitivity 
of our observable diminished significantly (see fig. 
||), but nonetheless remained sizeable. Enough to 
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Figure 2. Acoplanarity distribution of the p + p~ 
decay products' in the rest frame of the p + p~ pair. 
A cut on the differences of the tt^ to n° energies 
defined, in their respective replacement r^ rest 
frames, to be of the same sign is used. All detector 
smearings are included. The thick line denotes 
the case of the scalar Higgs boson and the thin 
line the pseudoscalar one. 
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guarantee the meaurement of Higgs boson parity 
to a confidence level greater than 95%, if typi- 
cal assumptions on linear collider luminosity, and 
(120 GeV) Higgs-boson production mechanism 
are made. Our method proved quite stable, with 
a degrading assumption on angular and energy 
resolution of the detector. Even a reduction by a 
factor of 2-3 with respect to typical assumptions 
used in linear collider designs did not destroy the 
method. To preserve its sensitivity, the angular 
resolution must be < 0(m p /E p ). 

3. MC-TESTER 

In the previous case the HEPEVT data structure 
was used to search for possible hard processes to 
calculate the spin state of the decaying r. One 
can search over an event record such as HEPEVT 
for other information as well. An essential step 
before attempting any sizeable rebuilding of the 
TAUOLA library is to automate at least some of 
its tests (such rebuilding may soon become nec- 
essary because of improving r-lepton data) . Such 
a test-program may be useful for other appli- 
cations, not necessarily related to generators of 
r-lepton decay, as well. That is why we have 
developed an individual program for that pur- 
pose. The idea behind MC-TESTER Q is quite 
simple: the user loads an extra library and, after 
the generation of every event with his generator, 
calls MC-TESTER, specifying the identifier for the 
particle to be searched for over the whole event 
record and studied. A data file from such a run 
is formed. 

In the second (analysis) step, data files from 
two runs of different Monte Carlo programs can 
be compared. The output is given in the form of 
a I^Tj^X file, which includes a table of all decay 
modes found(see fig. ||). The table is not only 
a list of these modes, but also, of the branch- 
ing ratios calculated from the two runs. Number 
quantifying maximum of shape differences is also 
given. This number is calculated as a maximum 
of the shape difference parameters (with a part of 
the code easy to identify) for all invariant mass 
distributions of the decay channel under study. 
We choose to take into consideration invariant 
mass of every subset of the decay products for 
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Figure 3. Second page of MC-TESTER booklet pro- 
duced at analysis step. 

the decay channel. The table is followed by a 
booklet of these invariant mass distributions (see 
fig. ^ for an example of an individual plot), which 
are grouped into separate chapters for every de- 
cay mode. 

4. TAUOLA 

As was already discussed in Victoria , many 
options of physics initialization are now available 
for TAUOLA (they can be constructed from the files 
included in the program distribution package). 
They often differ by the models used for fits to 
the experimental data, and/or by the data them- 
selves. Let me stress that, in general, the best 
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(1) 



returns 0.0271364 with the numerical coefficient 
Cn = 1261/120 — 7r 2 (for compatibility, we 
use as in f| m T = 1.7842, /„ = 0.0933 and 
cos^c = 0.975). A later calculation (formula 33 
in @) finds C N = 1009/96 - tt 2 and result is 
0.0272249. With the computer power available at 
the time of j^], we could not distinguish between 
the two options. At present, such an exercise is 
straightforward and we find that TAUOLA returns 
0.0272060 ± 0.000011. We can see, that the value 
of Cjv from ref. |23|] is statistically favoured over 
the one from ref. l20| . 



Figure 4. Example of individual plot from analy- 
sis step of MC-TESTER. 

parametrizations will always be in the hands of 
the leading r-lepton experiments of the time. Our 
original version of TAUOLA initialization was in- 
deed meant to be used only after re-initialization 
of its physics content by the users. Nonetheless 
for experiments in the design phase or for which 
r-leptons physical processes help in other stud- 
ies, such a solution is not good. We are gradually 
including new options into TAUOLA. 

In the last two years, parametrizations of the 
hadronic form factors |L7j based on the Novosi- 
birsk low energy e + e~ annihilation data ]l8| were 
introduced into TAUOLA as a possible new option 
for r~ — > 7r~ ir~ ir + ir a 'v and r~ — > 7r _ 7r 7r 7r i/ 
decay modes |l9| . This is the only published con- 
tribution documenting upgrades on the physics 
content of TAUOLA over the last two years. 

4.1. Comment on one of the old tests 

In ref. || (formulae (37) and (38)) we have 
reported that, for r decay into 47r, when the 
chiral limit is taken (the case we used for some 
tests of phase-space generation), the results of 
TAUOLA agree with those of ref. |20 , only up 
to an overall factor^] |^ cos 2 9 C . In the case of 
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jblems with normalizatio n ( or/and 
J were again addressed J21|. The 



5. SUMMARY 

Let me recall finally the main developments 
of the TAUOLA package in the last two years: 
as stand-alone generator of r-lepton decays, 
and for its interfaces to other programs. The 
progress on r-lepton decay alone consists of a new 
parametrization of form factors for the 4-7T decay 
channels. This parametrization, mainly based on 
the Novosibirsk data, form a step toward TAUOLA 
as a framework, where different data and mod- 
els can be compared. In such multi parametriza- 
tion approach the appropriate tools are necessary. 
The MC-TESTER, helpful in different types of works 
for programs such as TAUOLA became available. It 
can also be useful, for example, to check if the 
interfacing of TAUOLA to other programs did not 
corrupt its initialization. The family of TAUOLA 
interfaces was enriched with the new option for 
its universal interface. It now provides complete 
spin correlations for Higgs boson decay into pair 
of r leptons. This option was instrumental in 

authors of that paper worry that it may be of grave conse- 
quences for TAUOLA as well. In practice, even if the normal- 
ization of the chiral current was wrong (by the same factor 
in the Monte Carlo and in the testing function), it would 
affect neither technical tests nor actual results of TAUOLA 
generation, since a different current is used then. Note 
that in cases such as CLEO initialization fea, the current 
is normalized directly from the r data collected by the ex- 
periments; the normalization originating from the model 
disappears completely. This is why the issue is only of 
marginal interest for TAUOLA as it is now. 
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designing an observable that is potentially useful 
to measure Higgs-boson parity at future experi- 
ments. 

Let me finally acknowledge people and collab- 
orations who contributed to the present shape 
of TAUOLA and related projects: A. Bondar, S. 
Eidelman, P. Golonka, S. Jadach, M. Jezabek, 
J.H. Kuhn, A. Milstein, T. Pierzchala, E. Richter- 
Was, N. Root, B.F.L. Ward, M. Worek, and 
ALEPH, CLEO (in particular A. Weinstein), Del- 
phi, Opal, L3 collaborations, the Karlsruhe the- 
ory group, and others. I would like to thank A. 
Pich for discussion. 
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